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ABSTRACT 

We report the discovery of a variable object in the Hubble Deep Field North (HDF-N) 
which has brightened, during the 8.5 days sampled by the data, by more than 0.9 mag 
in Igi4 and about 0.7 mag in V606> remaining stable in B450. Subsequent observations 
of the HDF-N show that two years later this object has dimmed back to about its 
original brightness in Igi4. The colors of this object, its brightness, its time behavior 
in the various filters and the evolution of its morphology are consistent with being a 
Type lb supernova in a faint galaxy at z ^0.95. 
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1 INTRODUCTION 

The discovery of supernovae (SNe) in the early universe is 
of great interest because they can provide a wealth of infor- 
mation about cosmological parameters and the cosmic star 
formation history. It is now believed that the star forma- 
tion activity in the universe very likely started in small ob- 
jects which later on merged to form larger units (Couchman 
& Rees 1986, Ciardi & Ferrara 1997, Haiman et al. 1997, 
Tegmark et al. 1997, Ferrara 1998). Unless the IMF at high-z 
is drastically different from the local one, some of the formed 
stars will end their lives as SNe. Detecting high-z SNe would 
be of primary importance to clarify how reionization and re - 
hating of the universe proceeded (Siardi & Ferrara 1997), 
and, in general, to derive the star formation history of the 
universe (Sadat et al. 1998, Madau, Delia Valle & Panagia 
1998) and pose constraints on the IMF and chemical enrich- 
ment of the universe. 



Gre a t effort has been spent in this search ( Garnavich ei 
1998 ; Perlmutter et al. 199S ) and many SNe have been 



found up to a redshift of z = 1.20 ( Aldering 1998) when the 
unive rse had only about ha lf of its present age. The HDF im- 
ages ( [Williams et al. 1996 ) are the among the deepest taken 
to date and in principle could contain SNe up to z ~3. Two 
SNeQwere actually detected by Gilliland & Phillips (1998) 
and Gilliland, Nugent & Phillips (1999) by comparing the 
primary HDF-N data with second-epoch images taken two 
years after. The primary observations of this field were taken 
using four optical filters centered at 3000A (U300), 4500A 
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(B450), 6060A (V 606) and 8140A (I814) during a time span 
of about 10 days in December 1995. Although the distribu- 
tion of the data along this period is not uniform, the overall 
time coverage is good enough to detect objects with signif- 
icant variations on time scales of a few days. High redshift 
SNe cannot be identified near their maximum in such a short 
time span because they vary too slowly, but soon after the 
explosion they evolve fast enough to be detected. 

How many SNe can be expected in the HDF-N? This 
number can be estimated using the computations by Marri 
& Ferrara (1998) and Marri et al. (1998) for future Next 
Generation Space Telescope (NGST) surveys. Scaling their 
results for a flat CDM+A with f2 M = 0.4 to the HDF-N area 
(~ 5 sq. arcmin) and assuming a surveying time of ~ 8 days 
(see below), we obtain an expected number ~ 0.34 SNe from 
massive stars (SN types lb /II) . As a comparison, scaling the 
analytical estimates in Miralda-Escude & Rees (1997) to the 
HDF-N, one obtains a similar expected rate of 0.17 SNe in 
the HDF-N. These values are large enough to encourage a 
new analysis of the HDF-N. 



2 OBJECT SELECTION, PHOTOMETRY AND 
MORPHOLOGY 

The original observations of the HDF-N in each filter con- 
sisted of about 300 images taken in several (from 9 to 11, 
depending on the filter) positions on the sky (dither posi- 
tions). We have divided these images into a few subsequent 
groups, three for V606 and Isi4 and two for B450. The U300 
band, which is intrinsically less efficient, was not considered. 
The HDF-N was observed again two years later in U300 and 
Igi4 and the latter image (14 in Table 1) is deep enough to 
be used for this project. 

Where possible, images taken in the same position on 
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the sky were not split into different groups so that we could 
make effective use of those combined by the HDF working 
team (Williams et al. 1996) for each dither position, and 
made available on the web. Some dither positions in V606 
and i8i4 contain images too widely separated in time and 
therefore have been split. In this case the frames were re- 
duced by the automatic pipeline provided by the ST-ECF 
web site. In all cases warm and bad pixels were rejected 
using the standard routines. 

Given the small number of images, between 3 and 6, 
present in some of the fin al groups, it was not possible to 
use the drizzle algorithm (Fruchter et al. 1997) often used 



for the WFPC2 data reduction. We have chosen to resample 
all the images to a pixel of 0.05 arcsec (half of the original 
WF pixel scale) and use the IRAF task IMCOMBINE to 
make the final combinations while rejecting deviant pixels. 
This procedure is more efficient than drizzle in rejecting any 
residual warm and bad pixel; this is especially true for the 
final V606 and Igi4 images which contain 5 dither images 
each. Table 1 lists the properties of the resulting images 
with their time coverage and limit magnitude. 

Two tests were performed to check the data reduction 
results: i) the flux of a few compact sources from the list 
in Mendez et al. (1996) was measured in each of the result- 
ing images to verify the constancy of the photometry; ii) as 
discussed, for example, in Ferguson (1998), the effective lim- 
iting magnitude depends on the object size as much as on 
its total magnitude. The detection limits for point sources 
in the images were measured by adding stars and detecting 
them using FOCAS (Valdes 1982): the derived values for the 



80% completeness level are in good agreement (within 0.1 
mag) with the values in Williams et al. (1996) once their 
10(J limits in an aperture of 0.5 arcsec are scaled for the 
exposure time and to about 3<r in 0.30 arcsec, an aperture 
yielding the highest SNR for point source photometry (eg., 
Thompson 1995). These values are listed in Table 1; all mag- 
nitudes are in the AB system. 

The three combined V606 images (VI, V2 and V3), 
showing the best sensitivity (3a limits between 29.1 and 
29.7) and time coverage (about 8.5 days), were examined for 
variable objects having a monotonic trend, either brighten- 
ing or dimming. A few interesting objects were selected; the 
most remarkable one is in chip 2 of the WF camera and is 
present in the Williams et al. (1996) catalog with the entry 
number 584.2. The J2000 coordinates of this object (2-584.2) 
are 12:36:49.357 +62:14:37.50. In the total images, this ob- 
ject has flat B450, V606 and Igi4 colors and is undetected 
in U300; however it cannot be classified as a U300 drop-out 
being too faint to show a strong enough break between U300 
and B 4 5o- 

Table 1 and Fig. 1 show the time evolution of the pho- 
tometry of 2-584.2. Its magnitudes in the various images 
were measured inside a circular aperture of 0.3 arcsec, cor- 
responding to about twice the PSF FWHM, and corrected 
to an infinite aperture. This object shows a strong brighten- 
ing both in V606 and Igi4, while the photometry in the B450 
band is consistent with a constant value. The errors shown 
in Fig. 2 are derived under the assumption of Poisson noise 
inside the photometric aperture. 

We have computed the statistical significance of the de- 
tected variation. For the V606 band we consider the two 
differences V2-V1 and V3-V2; in the Isi4 band, as the ob- 
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Figure 1. Photometry of 2-584.2 in B450, V606> and Igi4 as 
a function of the (observer) time; day=0 corresponds to the first 
observation. Each point shows the 1 <r photometric errors and the 
time span covered by the data; the dotted lines are the 3 a limits of 
each image. The the arrows in the Igi4 panels show the 3 <r upper 
limits. The solid thick lines are a fit obtained from a K-corrected, 
time-dilated SNIb light curve at z = 0.95 with t maa; =34.0, E(B- 
V)=0.04 and £7o=0.1. The spectrum and the time evolution of a 
SNIb simultaneously fit the magnitudes of the variable object and 
its evolution in the three filters. The underlying galaxy results to 
have B45o=28.70, V606=29.10 and Igi4=29.15. Also shown are 
the curves for a SNIa at 2 ~ 1.3 (short dashed) and a SNII at 
z ~ 1.1 (long dashed) which do not provide an acceptable match. 



ject shows no significant variation among II, 12 and 14, we 
coadd these three images and compare the result with 13. 
By comparing the differences in the photometry with the 
quadratic sum of the errors, we find a joint probability for 
these three variations to arise from noise of 1.2x10 . In 
the HDF-N t here are about 1200 objects between Igi4=27 
and Igi4=29 (Williams et al. 199C); therefore we expect one 
such spurious event only every 70 Hubble Deep Fields. 

Fig. 2 shows the images of this object in the various 
bands [] and its brightening in Vsoe and Igi4- In the total 
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Table 1. Grouping of the images and properties of 2-584.2 



Filter 


Image 


Dither posj^j 


Exp. time 


Time range 


Limitf] 


2-584.2 








(hours) 


(days) 




Mag 


snr| 


1450W 


Total 




33.5 


0.00-9.21 


29.72 


28.73 


5 1 




Bl 


6A+9+4+5+3+8+1+2 


28.2 


0.00-3.15 


29.64 


9R 7^ 


7 r i 




B2 


7+6B 


5.3 


8.34-9.21 


28.72 


28.41 


3 8 


F606W 


Total 




30.3 


1.48-10.07 


30.07 


98 78 


7 9 




VI 


3+2+6A+4A 


5.6 


1.48-5.56 


29.16 


29.14 


3.1 




V2 


9A+5A+7+1+11+10 


14.7 


6.43-8.51 


29.68 


28.93 


6.0 




V3 


8+6B+4B+9B+5B 


10.0 


8.57-10.07 


29.47 


28.53 


7.2 


F814W 


Total 




34.3 


3.21-10.10 


29.46 


28.63 


4.8 




11 


2+1+6 


10.0 


3.21-4.35 


28.80 


>28.80 


2.8 




12 


4A+9+3+5A 


14.5 


4.69-6.30 


29.00 


28.76 


3.8 




13 


8+7+11+4B+5B 


9.8 


6.63-10.10 


28.79 


28.05 


6.0 




14 




15.6 


735-738 


29.00 


>29.00 


2.9 



a Dither positions as in Williams et al. 1996. When present, A and B refer 

to the first and second part of the split positions 

b 3<t limit for point sources 

c Signal-to-noise ratio of the photometry 



images the object is marginally resolved in B450 and V606, 
whereas in Isi4 it is consistent with a point source. The im- 
ages also give the impression of a sharpening of 2-584.2 with 
time: the object seems more extended at the beginning of 
the observations than at the end, as if a bright core were 
emerging in V606 and Isi4. Its faintness prevents us from 
studying its morphology in detail; nevertheless, its exten- 
sion can be roughly measured by fitting circular gaussians 
to each image. The results of this procedure supports the vi- 
sual impression of a sharpening of the object with time: the 
FWHM of the best-fitting gaussian goes from 0'.'33±0'.'08 
to 0'.'18±0'.'02 in V 6 oe and from 0'.'40±0'.'10 to 0'.'19±0'.'02 
in I814 while the value expected for point sources is about 
0'.'15. 



3 THE LIGHT CURVE 

Among the objects that could show variability at this high 
galacti c latitude and faint flux level. AGNs and SNe are 



the mo >t plausible candidates. The shape of detected time 
variation and the possible development of a bright core in 
an underlying galaxy strongly suggest the identification of 
2-584.2 with a SN observed soon after its explosion. 

The expected evolution of the apparent magnitude of a 
SN can be derived from template light curves, spectra and 
absolute magnitude (as discussed below) once its distance is 
known and K-corrections (due to the narrowing of the filter 
passband in the restframe of the source and redshifting of 
the emitted photons from the source to the observer) are 
computed. For high redshift objects the time in the observer 
frame, t bs , is related to the time in the supernova rest frame 
by tots — (1 + z)t re3 t; this produces a time dilation of the 
high redshift SN light curve. 

Once the SN template is given, four parameters are left 
to simultaneously fit the data in the various filters: the SN 
redshift, z, the time of the maximum light, t max , the color 
excess due to dust, E(B-V), and the cosmological matter 



density parameter, Slo (we assume A = 0) . Some uncertainty 
remains associated with this fit as pre-maximum light curves 
and spectra (particularly in the UV) are generally not very 
well determined because they are available only for a handful 
of objects. The SN flux must then be added to the host 
galaxy flux in each filter. The observed time variation of 2- 
584.2 shows that the SN becomes quickly dominant in V606 
and I814 while the galaxy produces most of the B450 flux. 

We now consider in turn the possibility that the de- 
tected variable source is a Type la (SNIa), a Type II (SNII) 
or a Type lb (SNIb) supernova. SNIa have been shown to 
be reliable standard candles, as their intrinsic luminosities 



can be accurately determined (Phillips 1993). Pre-maxi mum 



optical (CTIO ) and UV (IUE) spectra of SN1990N (Lei- 



bundgut 1991), show that the flux drops off sharply below 
2600A; therefore, as soon as z Si 0.7, their B450 flux in neg- 
ligible as the cut off is redshifted into this band, while the 
mere detection of variation in the V606 band puts the up- 
per limit z ^ 1.6. The lumin osity of a SNIa evolves in time 



according to a light curve (Leibundgut 1988; Doggett & 



Branch 198E ) which shows a fast rise to the maximum (3.6 
mag in 15 days), with slight differences between the (rest 
frame) B- and V-band. When fitting the data for 2-584.2 
with the time-dilated, K-corrected light curves, we find im- 
possible to simultaneously match the Veo6 and Isi4 data 
points (in this case the B450 data only give information on 
the galaxy magnitude), as shown in Fig. 1. Since the object is 
relatively faint, this implies high SNIa redshifts (z ^ 1.3) for 
which colors cannot be reproduced by using appropriate K- 
corrections. We conclude that identification of 2-584.2 with 
a SNIa can be safely ruled out. 

We repeat the same procedure for SNII. These objects 
cannot be used as standard candles as their peak abolute lu- 
minosities are kn own to cover a wide range, from Mb ~ — 14 
to M B ~ -19 dPatat et al. 1994|). SNII are usually di- 



vided into two classes (Doggett & Branch 1985), one (SNII- 
P) showing a slow pre-maximum brightening and a plateau 
in the after-maximum decline, the other (SNII-L) a faster 
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Figure 2. Images of 2-584.2 in the B450, V606 and L314 bands. The three images on the left are the total HDF-N images by Williams et 
al. 1996 with a pixel scale of 0.04 arcsec/pixel. On the right the time sequence in each filter with pixel scale of 0.05 arcsec/pixel showing 
the brightening of the object in V606 an d Isi4 and its following dimming. Each image is labelled by the filter and the time range of the 
observations in days. The size of the central squares is 1 arcsec. 



brightening and a linear decline. In both cases the pre- 
maximum spectra can be approxima ted by a black -body 
with a temperature T B b ~ 25000 K (|Kirshner 199(i|) with- 



out any UV cut-off. By comparing the expected light curves 
with the data (see Fig. 1) we can safely exclude both classes 
because (i) their brightening is too slow and (ii) their blue 
spectrum makes them too luminous in B450. 

SNIb yield a good agreement with the data. These ob- 
jects closely resemble SNIa in terms of time evolution and 



spectra but are typically 1 -2 mag fainter (Wheeler & Lev 
reault ] 985| ; Kirshner 199C) and show a larger spread in the 



maximum brightness. While SNIa are found in all types of 
galaxies and derive from old stars, the SNIb are only de- 
tected near regions of active star formation and their pro- 
genitors should be young massive stars. As shown in Figure 
1, both the time evolution of 2-584.2 increasing from B450 to 
V606 to I814 and its apparent magnitudes are easily repro- 



duced by SNIb light curves ( Kirshner 1990] ): acceptable fits 
can be obtained for 0.90 < z < 1.02 which makes this ob- 
ject one of the most distant SN observed to date, while the 
best agreement is found for 2 = 0.95, tmaa;=34.0 ± 1 days 
(~ 12 rest-frame days after the end of the observations), 
l ow Qq value s and moderate reddening <E(B-V)< 0.05 
flSeaton 197&] ). This values of the reddening, corresponding 
up to about 0.29 mag of extinction in Veoe and 0.22 in Isi4, 
could also be accounted for by a SNIb fainter than the av- 
erage and without extinction. Lower values for the fitting 
redshift tend to select low values of f2o and low extinction, 
viceversa for higher redshifts. The galaxy, showing flat B450- 
V606 and V606-I814 c olors and luminosi tes indicating a SFR 
of about 0.2 Mq/yt (Madau et al. 1998), is consistent with a 
star forming dwarf. About 20% of the flux in the Isr4 image 
taken two year later (14) is still due to the SN. 

We therefore conclude that: (i) SNIa are too bright and 
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red and SNII are too slow and blue to be viable candidates 
for 2-584.2; (ii) a SNIb naturally reproduces the time evolu- 
tion, brightness and colors of the variable object. This type 
of study, when applied to high sensitivity, long time span 
future observations, could produce a significant sample of 
high-z SNe that can be succesfully used to constrain the 
star formation history and the geometry of our universe. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science IATfjX style file. 
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